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ABSTRACT 

The motion of compact objects in Galactic halo is investigated in the approximation of point 
mass under the assumption that the external halo is not gravitationally dominating and consists 
of test bodies. Then, the phase space distribution function is that of a spherically symmetric 
collection of confocal elliptical orbits of various energies and eccentricities, called Keplerian 
ensemble, and can be determined in a systematic way from the radial velocity dispersion pro- 
file of the halo objects. In particular, the anisotropy of the velocity dispersion tensor becomes 
a function of Galactocentric distance. We find that the observed motion of halo objects can be 
made consistent with Galaxy mass of 2.1 x lO^^M©. 



1 INTRODUCTION 



Galaxy is customarily modeled as a composition of several sub- 
systems described by simple mass profiles: the central bulge and a 
spheroid of stars and gas, the galactic disk and an extended exter- 
nal halo. The halo is assumed to consist of nonbaryonic dark mat- 
ter (NDM) mainly, whereas ordinary matter makes up only a small 
fraction of the halo. In this picture, NDM dominates the gravita- 
tional potential at larger radii. 

There are large discrepancies in the estimates of Galaxy mass 
and uncertainties in its distribution. The available data cannot dis- 
criminate between various mass models. Tenths of Galaxy mass 
estimates are given to chose from, depending on the particular halo 
model, its free parameters, and the choice of satellites assumed 
gravitationally bound to Galaxy. A good deal of this indeterminacy 
lies in the mass-anisotropy degeneracy characteristic of halo mod- 
els. This uncertainty concerns mainly the non-baryonic dark mat- 
ter (NDM) halo. There are opinions that the mass distribution is 
the most poorly known and c ontroversial Galactic characteristic at 
large distances (iBrown et~aLL201Q) and th at the halo issue is en- 
tirely speculative dPehnen & Binnevl[l998h . The problem of NDM 
concerns various scales of distances. In the Galaxy interior its pres- 
ence can be inferred from the rotation curve, the vertical gradient 
of the rotation and from gravitational microlensing (there are also 
attempts at direct detection of NDM in Earth based experiments), 
whereas at larger distances, NDM presence is inferred from the mo- 
tions of compact objects in the external halo. 

Except for the habitual argumentation for an extended and 
massive NDM halo, there is also a premise for not considering it, 
at least in the Galaxy interior. The high vertical gradient of ro- 
tation above the disk midplane can be e xplained in the thin disk 
model approximation ( Jalocha et al. 2010) (whereas spherical sym- 
metry spoils this agreement). Moreover, in a related finite width 
disk model, the amount of mass ascertained through microlensing 
inside the Sun's orbit is found consistent with the dynamical mass 
inferred from rotation in the same region dSikora et allboilh . The 
two observations suggest that the overall mass distribution might be 



flattened rather than spheroidal. Furthermore, the influence of inter- 
stellar magnetic fields on rotation of gas can be responsible for flat- 
tening of rotation curves a t larger radii beyond the visible galactic 
disks ( Battaner et al.l ll992), which is ignored in modeling of rota- 
tion curves. In this context it is worth to mention that the character- 
istic rise in the mass-to-light ratio profile (interpreted as a signature 
of NDM) co uld be reduced whe n the magnetic effect was taken 
into account dJalocha etal.]l201lh . What is more, recent indepen- 
dent experiments aimed at the search for dark matter consisting of 
weakl y interacting massive partic les, detect no trace of dark matter 
signal (lAckermann & et al.ll201ll) and exclude generic WIMP can - 
didates for dark matter (iGeringer- Sameth & Ko ushiapDasI l201lh . 
An experiment aimed at direct search of NDM also fails to find 
it ( Aprile & at aD[201lh . The null results of the experiments could 
be reconciled with the amounts of dark matter ascertained through 
gravitational lensing at the scale of galaxy clusters. It is possible 
that the clustering scale of dark matter is larger, in which case dark 
matter would tend to clump around large clusters of galaxies rather 
than around separate galaxies. 



Point mass approximation as a means to investigate motions in 
Galactic halo 

The asymptotic value of mass function of the Galactic gravitational 
potential (total Galaxy mass) can be obtained by investigating the 
motions of the most distant compact objects scattered in Galactic 
halo. They can be treated as test bodies moving in the background 
gravitational field. A primary quantity used to describe their motion 
in a statistical manner is the phase space distribution function. The 
radial velocity dispersion of halo objects and (normalized) number 
density are secondary quantities defined as appropriate integrals in- 
volving the phase space distribution function. 

Mass density of matter producing gravitational field is not nec- 
essarily directly related to the radial velocity dispersion of test ob- 
jects moving in that field. As an example of a system realized in 
nature for which this is true can be given the structure of Saturn 
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rings. The countless particles (of which the rings consist of) move 
as test bodies in the gravitational field of the planet and the rings 
do not contribute significantly to this field - the gravitational po- 
tential is unrelated to the density function of the rings in the Jeans 
equations. 

It cannot be excluded the possibility that the heavy NDM halo 
is not present in the Galaxy and this hypothesis should be tested. In 
this case, at radii sufficiently large, the gravitational field would be 
approximately that of a point mass - the contribution to the grav- 
itational field from higher multi-poles of the compact mass distri- 
bution (consisting of Galactic disk and a spheroidal component of 
which part is the central bulge and a halo of gas and stars) would 
only slightly affect orbits of external objects. These objects would 
not contribute significantly to that field and could be regarded as 
test bodies for which the Kepler problem applies with the central 
mass representing total Galaxy mass. Thus, when gravitationally 
bound to Galaxy, every compact object of the halo could be as- 
sumed to move in the asymptotic field as a test body on an elliptic 
orbit. One can study a collection of such orbits and calculate var- 
ious expectation values over concentric spherical shells centered 
at the origin. It is evident that the radial velocity dispersion of the 
objects or the velocity dispersion anisotropy would not be directly 
related to the gravitating mass distribution - owing to the small 
number of halo objects they could be added or removed without 
changing the field, whereas the radial velocity dispersion profile 
could be changed significantly. 

The idea of studying the motion of halo objects in the point 
mass approximation should not a stonish, as it is not new. It was 
considered already 30 years ago bv lBahcall & Tremainel (ll98l') and 
applied to several external galaxies. They proposed some estima- 
tors of galactic mass, all being an average of the quantity v\RIG 
over external halo objects, times a constant factor. The factor de- 
pended on the a ssumption about the mean square of eccentricity. 
A recent paper ('Watki ns et all 12010') offers in essence the same 
method as originally proposed bylBahcall & Tremaine (1981). The 
difference is the presence of arbitrary power of distance in the mass 
estimators (owing to the assumptions about the spatial distribution 
of halo objects and the velocity dispersion anisotropy). Using esti- 
mators of this kind is tantamount to assuming a very special family 
of phase space distribution functions, which is a strong constraint. 
Interestingly, they found that the Galaxy mass could be as low as 
4 X 10^^ M© (but, among tenths of other masses, they give also 
2.7 X 10^^ M© as their preferred result if Leo I is assumed gravita- 
tionally bound to Galaxy). 

Our approach to investigating the motion of halo objects in 
the point mass approximation is essentially different - it does in- 
volve estimators of this kind. We allow for arbitrary phase space 
distribution function, it is not a priori assumed. It is quite a generic 
function of two variables (energy end eccentricity) with the only as- 
sumption that the system is spherically symmetric. Expectation val- 
ues of various characteristic quantities such as the dispersion tensor 
anisotropy, the profile of the radial velocity dispersion, all become 
functions of Galactocentric distance. In particular, we abandon the 
constraining assumption of setting the anisotropy parameter as a 
global constant. Our model allows for the possibility that on aver- 
age the external orbits could be more elongated than internal ones. 
One of advantages is to avoid overestimation of Galaxy mass. We 
apply a minimizing scheme that enables to find the optimal form of 
the phase space function from the measurements of halo objects. 

The result of the present work is to show that the observed 
motion of halo compact objects (like stars, satellite globular clus- 
ters or dwarf galaxies) can be made consistent with a low Galaxy 



mass of about 2.1 x lO^^M© at a distance 180 kpc, much lower 
than predicted by other models, and only 7/5 times greater than 
that predicted in the disk model framework based on the Galaxy 
rotation curve (extending out to 20 kpc). With this mass Leo I and 
and Hercules cannot be gravitationally bound to Galaxy. 

The paper is organized as follows. In the next section we give 
simple estimates for Galaxy mass. This will be similar to what other 
papers offer in the point mass approximation. The third section con- 
tains our general solution to Jeans equation in a point mass field 
for a general spherically symmetric collection of confocal orbits of 
test bodies moving in the space limited by two concentric spherical 
shells centered at the origin. This requires considering a finite sup- 
port for the cutoff phase space distribution function and a related 
problem of choosing the appropriate set of independent variables 
on a fragment of a the plane of energy and eccentricity. Also some 
expectation values over thin spherical shells are defined which are 
indispensable in determining the mean values of various observ- 
ables as functions of the galactocentric distance. In the fourth sec- 
tion we determine the radial velocity dispersion for compact objects 
in Galaxy halo, then we discuses briefly a minimizing scheme for 
finding the phase space distribution function in the framework of 
Keplerian ensemble, and finally we apply our method to find the 
Galaxy mass in the point mass approximation. Then we give con- 
clusions. 

Throughout this paper we assume the following parameters: 
Rs = 7.62 zb 0.3 2 kpc for Sun's distance from the Galactic cen- 
ter (Eisenh auer &at alll2005h. V. = 208 ± 16km/s for the local 
disk rotation speed jBrainerdl [2007h and {U,V,W} = {9.96 di 
0.33, 5.25 ± 0.54, 7.07 ± 0.34} km/s for the components of the 
velocity vector of the Sun with respect to the local standard of rest 
( Aumer & Binnev 2009) (for a summary of other measurements, 
see (iFrancis & Andersonll2009l) ). 



2 SIMPLE ESTIMATES OF GALAXY MASS 

In the case of spherical systems one can look for a relation M(r) = 
G~^r a{r){vr{r)) between mass function M(r) and radial ve- 
locity dispersion function {vl(r)) determined over thin spherical 
shells. The coefficient a{r) is a dimensionless function of galac- 
tocentric distance r. The connection between M(r) and {vr{r)) is 
model-dependent and the predicted asymptotic mass can be differ- 
ent for the same halo measurements. 

The assumed mass function in NDM halo models is related 
to the theoretical {vr{r)) profile through Jeans and Poisson equa- 
tions. There is no freedom in modifying the phase space distribu- 
tion function (modulo free parameters) since it is assumed. The 
free parameters are found by minimizing a measure of distance be- 
tween {vr(r)) and halo measurements. The parameters, and so the 
total mass, are not uniquely determined, the same dispersion pro- 
files could be explained by models with different total masses (this 
is the mass-anisotropy degeneracy). In the case of a point mass 
source, the coupling function a{r) is a functional of eccentricity 
and energy over various orbits of test bodies, whereas M(r) is 
constant. There are possible infinitely many functions {vr{r)) in 
the same gravitational field. This freedom would be limited by im- 
posing an ansatz for the number density of compact halo objects, or 
by assuming the dispersion anisotropy parameter as constant, etc. 
Starting from section [3] onwards, we will refrain from making such 
constraining assumptions. 
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2.1 Galaxy mass from escape velocity 

At distances from the Galactic center sufficiently large, the gravita- 
tional field can be approximated by that of a point mass M located 
at the center and regarded as the total gravitating mass of Galaxy. 
For a distant test body bound in such a field 
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e2 < 1. 



2GM 

Here, ^ e < 1 is the eccentricity of an elliptical orbit, E and J 
are the energy and angular momentum (per unit mass), and 

2^J^ 2 . ^ G^M^ 



1 + 



> ^ ^ 



Hence, there is an upper bound for the absolute value ofVr = f 
possible for a remote halo object, Vr^ < provided the object 

is gravitationally bound. One can try to estimate the lower bound 
for the Galactic mass by inverting the inequality. This gives us the 
escape mass observable in spherical symmetry 



M{r) 



2G 



< M. 



A histogram showing the number of catalogued halo objects with 
A4 not greater than some limit is presented in figure [T] 
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Figur e 2. Galaxy rotation curve. A unified rotation curve from dSofue et 
"2009*) prepared at Rs = 200 kpc a nd Vs = 200 km/s [ circles+error 
bars] and their model rotation curve ( iHonma & Sofue|[l996h [solid thick 
line]. The extended rotation curve in disk model [thin solid line] 




Figure 3. Keplerian mass function rv'^/G for Galaxy rotation data [circles] 
and a linear fit M(r) = 9.72 x 10^ R [solid line]. To compare with, the 
disk model mass function has been also shown [dotted line]. 



Figure 1. Escape mass histogram for halo objects. It shows the number of 

r (r) 

objects with — not greater than a given value (in solar units). 



Among halo objects there can be some that are not gravita- 
tionally bound to Galaxy and should be excluded from the analysis. 
Leo I with the highest M = 10.3 x 10^^ M© is such an object (this 
will be justified later). As so, the estimate for the Galaxy mass is 
3 X 10^^ M© at R ^ 130 kpc, provided that Hercules is bound 
to Galaxy on almost radial orbit. This value is to be compared 
with the integrated Galaxy dynamical mass of 1.44 x lO^^M© a t 
R ^ 20kpc predicted in the global disk model (Ij alocha et al.l201Qh 
for the rotation curve shown in Fig.|2] This value is smaller than that 
for a spherically symmetric mass model in Fig.[3] which gives the 
dynamical mass estimate of about 2.0 x 10^^ M© at R ^ 20 kpc. 
The lower bound for total Galaxy mass at R ^ 20 kpc can be 
something in between 2 x 10^^ M© and 3 x 10^^ M©. Then only 
^ 30% of the total Galaxy mass would reside outside = 20 kpc. 



2.2 A simple mass estimate from Jeans theory 

For a spherically symmetric system in a steady state with vanishing 
mean velocities 



rp,r 



+ 2/3 = 



r 



from Jeans equations. On assuming that (vr) is almost constant at 
radii sufficiently large, that the transverse and radial dispersions are 
proportional to each other so that p is constant, that r $ ^ —GM 
and p ^ po {ro/r)^~^^ (with small £ > to ensure integrability at 
infinity), then ar {vl) ^ GM if a = 4 + e - 2/3. We see that M 
can be made arbitrarily large if the velocity ellipsoid is sufficiently 
flat (/3 <^ 0). A reasonable assumption for the simple model is a 
marginally integrable halo (s = 0^) with isotropic velocity disper- 
sion (/3 ^ 0), then a{r) ^ 4: 

4r 
'G 



M 



(Vr). 



By applying this "flat" model outside r = 30kpc to the (asymp- 
totically nonconstant) dispersion profiles we present later (see 
Fig-HB, one can obtain the following estimate for Galactic mass: 
5.4 X 10^^ M© if we discard Leo /; 4.4 x 10^^ M© if we discard Leo 
I and Hercules; and 4.0 x 10^^ M© if we discard Leo I, Hercules 
and Pal 14. 

Corresponding to the above model is the following mass es- 
timator M — ^ . TiVrl using a number of N outermost halo 
objects. A similar class of estimators in the point m ass approxima- 
tion was considered in teahcall & Tremainelll98lh . The result is 
shown in Fig.|4las function of the number of most distant objects, 
assuming that Leo I is not bound to Galaxy. The function has a dis- 
tinct minimum and a following maximum in which region the num- 
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Figure 4. The result of applying the estimator of total mass M.n = 

N 

"^i^rf to outermost compact objects in Galactic halo shown 

i=i 

as a function of A^. X = LogioN, Y = Mn/'^q- Leo I is excluded 
from the analysis. The most distant object is Leo T at r = 417 kpc with 
rv'^/G ^ 2.5 X IQiiM©. 



ber of objects grows from a dozen to a hundred without a signifi- 
cant change in the estimated mass. The maximum 4.6 x 10^^ M© 
can be taken as a safe mass estimate. As it is seen from the di- 
agram, the same value of mass could be obtained by taking into 
account only 41 most distant objects, with the one closest to the 
center located at ~ 46.1kpc. A similar analysis with the exclusion 
of Leo I and Hercules gives 4.4 x 10^^ M© with 37 objects outside 
46.7kpc, and with the exclusion of Leo I, Hercules and Pal 14 gives 
4.3 X 10^^ M© with 26 objects outside 50.9kpc. With such masses 
both Hercules and Pal 14 with escape mass 2.95 X lO^^M© and 
2.12 X lO^^M©, respectively, could be bound to Ga laxy. (To com- 
pare with, for the sample of halo objects studied in dWatkins et al.l 
|2010^ the mass estimate would be 5.5 x lO^^M© without Leo I 
which compares to our 4.6 x lO^^M0;or8.5 x 10^^ M© withL^o 
/ included - in either case Leo I could not be bound to Galaxy). 

The conclusion is that high Galaxy mass estimates come from 
large values of the theoretical coefficient a, depending on models. 
As seen from table [T] the coefficient could be increased by a factor 
of 1.25 by assuming phase space distribution function /(e, e) = 
e^/^ (with constant anisotropy (3 = —1). This way Galaxy mass 
could be increased without limits. Our goal is quite opposite, we 
aim to reduce a{r) to a minimum admissible by the observational 
data using a model maximally flexible with the minimum number 
of assumptions. This cannot be achieved using simple models with 
constant (3. 



2.3 A hindsight for mass observable 

The time average of rv^ over an elliptical orbit with eccentricity e 
is proportional to the central mass M: = M, similarly, 

Y^^{-^)t — M. These two averages are equal when = 2/3, 

2 

in which case M = 3(^)t. Analogously, for equal mean values 
of vz. and vl in the plane of motion, 1 = , — f , 

2 

hence = | and |(-^)t = M. With this hindsight, we choose 
an observable which best suits our task at hand. This must be some 
average of rv^. For a collection of orbits, the mean value of {rVr)t 
can be determined knowing the number of orbits with particular 
energy and eccentricity. In the next section we shall deal with aver- 
aging over concentric spherical surfaces. 



3 SPHERICALLY SYMMETRIC ENSEMBLE OF TEST 
BODIES IN THE FIELD OF A POINT MASS 

The motion of a test body in a spherically symmetric potential $(r) 
is planar in a plane through the origin. The plane is determined by 
two integrals of motion - a unit normal to the plane. The additional 
two integrals of motion are energy E and magnitude of angular 
momentum J (both per unit mass). In terms of angular velocity cj, 
cj^ = + 0^ sin^ 0, the two integrals read 

In the particular case of Newtonian potential 

$ = , 

r 

there is an additional integral of motion indicating a fixed direction 
in the plane of motion. 

The condition for all r > with J ^ gives 

1 + ^f^2 ^ 0- We are interested only in spatially bound orbits. 
Two turning points will be possible (or Vr = 0) when E < and 
1 + ^f^'} < 1. Hence, elliptical orbits are possible only for pairs 
{E, J) such that 0^1 + < 1. Elliptical motion can be thus 

uniquely determined by specifying five numbers: three Euler angles 
describing orientation of the orbit in space and a pair of numbers 
(e, e) defined by 

Here, R is some unit of length, e is the eccentricity and e is a mea- 
sure of energy determining the length of the large semiaxis, which 
is ^. The turning points are ^ (1 ± e). The dimensionless param- 
eters (e, e) play the central role in further considerations. 

To find various expectation values for a spherically symmet- 
ric collection of confocal ellipses, it suffices to know a distribution 
function describing the number of ellipses of a given eccentricity e 
and energy — e. It will be related to the distribution function /(r, v) 
in the /x-phase space. We assume that / can be expressed through 
first integrals e and e. Then, / is stationary and satisfies the nec- 
essary condition dtf + dqif -\- dqifd^^ = for a distribution of 
a collisionless system. This is the wording of the Jeans theorem in 
our situation. 

As an example, consider an almost axisymmetric collection of 
confocal ellipses randomly distributed in a plane (see Fig.O. One 
can assign to it a substitute smooth density function describing the 
expected number of bodies in a given surface element, assuming 
that every ellipse corresponds to a single body. In finding the den- 
sity one should take into account the fact that the probability of 
finding a body at a given point on its trajectory is proportional to 
the fraction of time spent in a given arc element. Given a family of 
orbits one can calculate expectation values for various observables 
at some distance from the center by taking into account the num- 
ber of various ellipses crossing a circle concentric with the origin, 
which all become functions of the circle's radius. The most impor- 
tant observables are: the radial velocity dispersion, the mean ec- 
centricity of orbits intersecting that circle, or the mean value of the 
velocity anisotropy parameter. Similarly, one can consider a spher- 
ically symmetric collection of confocal ellipses in space, which is 
the subject of the next section. 
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Figure 5. A random collection of confocal ellipses in a plane generated so 
as all ellipses are entirely contained within some annulus (here 1 < r < 12) 
with uniform probability for all admissible pairs (e, e). 



3.1 Expectation values for Keplerian ensemble 

The integrals of motion e and e can be treated as new indepen- 
dent phase space variables. By making a transformation from or- 
dinary spherical coordinates r, ^, 0, Vr,ve, to new coordinates 
0, 0, e, e, ip of the form 



Ru, 



GM f 2 _ 1 



2e 



^^^S (sin ^, COS 



the original volume element r^drAsin dOAd(l)A2-^f^AdveAdVci) 
(factor 2 to account for two signs of Vr) is transformed (up to a 
constant factor) to 

du A sin ^ A d(/) A dV^ A de A e de 



For we have already assumed spherical symmetry, the angles 0, 
and ip can be integrated out. The form of the volume element 
defines the integration region. 

In finding expectation values as function of radial distance, in- 
tegration must be carried out over all ellipses intersecting a spher- 
ical thin shell of some fixed radius r centered at the origin. Given 
< r < oo and ^ e < 1, we have ^ e 
orbits. Hence, we arrive at the distribution integral 



^ for such 

^ 2u 



where 



r 

/ f{f,v)d^fd^v oc 



ede / de 



f 

Jo 



du Uu [f] 



(l±e 
V 2u 



For physical reasons, it is convenient to assume that all orbits 
of the ensemble are contained entirely within a spherical annulus 
Rua < r < Rub (in between two limiting spheres of radius 
ra — Rua and = Rub). This requirement imposes restrictions 
on the integration region in the phase space, changing considerably 
the support of /(e, e). Then, Ua < and < Ub. Hence, 



l + e 
2nf, 



< e < 



l-e 

2Ua. 



and thus also ^ e < 



< 1. As a re- 



sult, the integration region has been shrunk to a quadrilateral abcB 
as shown in figure [6l However, integration requires splitting of the 























^^^^^^^^^ B 




_ c 



Figure 6. The points on the plane (e, e) represent confocal elliptical or- 
bits of eccentricity e and energy — e. The ABC triangle with vertices A 
(0, T^), B ( ""b-na 1 \ ^ (o, ^) comprises all orbits con- 



tained entirely between two spheres of radii Ua and Ub and centered at 

2u 



the focal point. The two families of lines e : 

l-e 
2u 



through point (-1,0) 



and e = V^/^ through point (1, 0) cross with each other forming a new 



coordinate system (o;, /3) on the plane (e,e). When a G (^^,1^ and 
/3 G (^1, the coordinates cover only a quadrilateral abcB with ver- 
tices a f ^^), b fo, :^), c f and B as required. 
Quadrilateral abcB is the locus of all orbits (e, e) which cross a sphere of 
some radius Ua < w < at least once. 

i^-dependent integration region into three parts. To overcome this 
difficulty we introduce a mapping (a, [3) (e(a, eu{a, [3)) 
from a rectilinear region with cartesian coordinates (a, /3) to a in- 
dependent irregular region abcB. This leads to the following coor- 
dinate change in the reduced phase space (e, e, in) in integral Vu [f] ' 



7^- < a < 1, 1< P < -r-, 0<ua<u<ui,. 

The interpretation and the origin of this coordinate change is de- 
picted in figure [6] We get finally. 



f{f,v)d^fd^v oc / du\/uiJ.u[f] 



where 



1 u/ u 

J VT^J 



dp 



13- 



1 



Now, given /(e, e), all expectation values can be in principle deter- 
mined. In calculating them the following expressions will be useful 

..2 o 



rVr 
GM 



GM 



a + /3 ' GM a + /3 + a 

For the purpose of further applications the most important are av- 
erages over concentric thin spherical shells. The mean value of a 
function g defined on a spherical shell of radius r is 



{9)r 



fJ^u[f] 



The above expression for {g)r is quite analogous to a conditional 
probability 'A provided that B'. Consequently, the average over all 
spherical shells can be calculated from 

/ d^ {g)ufiu[f] 



(9) = 



J du/Ju [f] 
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As a simple example, consider /(e, e) = Ce^^^ in the limit Ua 



0, Ub oo. Then, 



{rvf.)r 



= I , the number density must behave 



GM ~ 4' 

as r""^ (the number of particles diverges logarithmically). With the 
same function / but with different support so that all orbits with 
u < Ua = 1/24 or u > Ub = 24 were not present in the system, 
the density profile and averages change considerably (see Fig.[7]). 
See also table[T]and figure Fig.[8]for various other /'s. 



p 




Figure 7. Various observables averaged over spherical surfaces shown as 
functions of the spherical radius with phase space distribution function 
/(e, e) oc e^/^. The expectation value for the density of objects p is in- 
tegrable to unity. First column shows expectation values with the support 
1/12 < li < 24 and the second column with 1 < w < 24. 



The lesson we gain is that in the halo modeling we should 
assume / ^ e^^^^'^ase 0,^ > (for integrable number 
density). In addition, we can assume i^a > to eliminate all orbits 
with pericenter smaller than Ua . Physically, this condition allows 
for considering only bodies moving sufficiently far from the center. 
In conjunction with the previous condition, the number density will 
be globally integrable. The former condition can be released by 
considering orbits for which u < Ub always. Then, the cumulative 
number of objects will start from zero aiu = Ua and grow until it 
saturates Siiu = Ub. 



4 APPLICATION TO GALACTIC HALO 

We will be using averages {rVr)/G over (reasonably thin) concen- 
tric spherical shells as the basic observable. It is the simplest one 
with the dimension of mass formed out of pairs (r, Vr) measurable 
for halo objects. It resembles a size independent time average for a 
single closed orbit {rVr)t — \GMe^ . In order that {rvl) could be 
determined from measurements as a function of r with a satisfac- 
tory resolution, a reasonable number of compact objects in Galaxy 
halo is needed in every of a number of sufficiently thin concentric 
spherical shells centered at the origin. 

To prepare the {rVr)/G pro file we used position -v elocity 
data from the following papers: IClewlev et all (l2004l) : iMated 
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Table 1. Averages over an infinitesimally thin spherical shell of radius r 
for a spherically symmetric ensemble of confocal elliptical orbits described 
by phase space distribution function /(e, e). These values are only for il- 
lustration (the number density p of test bodies is nonintegrable). Parameter 
/3 is the anisotropy parameter of the velocity ellipsoid. With these simple 
/(e, e)'s all the quantities (except p) are independent of r. 



Wilhelm etaD (Il999|); IStarkenburg et all (l2009h : 



iDohm-Palmer et all (l200lh : iHanisI fmd) and recalculated to 
epoch J2000 when necessary. In addition, we included objects 
like Ursa Major I and //, Leo IV, Leo T, Coma Beren ices, Canes 
Venatici I and //, Hercules iSimon& Gehal l2007h. Bootes L 
Willman 1 (Martin et al. 20^, B ootes II (iKoch et al.ll2009l) . Leo V 
dBelokurov et al. 2008), Segue I dOeha et alJl2009h . and Segue II 
jBelokurov et al. 2009). In total, we used 1361 compact objects in 
our sample. 

The sample of halo objects is shown on the (r, rvl/G) plane 
in Fig.|9] It is not representative and is insufficient to determine a 
reliable profile of {rvl) at larger radii. Most of objects are located 
in the Sun vicinity along a few directions only and the statistics 
is very poor at the most crucial for Galaxy mass determination re- 
gion of the outermost radii where only several objects are present 
spread over a wide interval of Galactocentric radius, see Fig.fTOl 
The problem of lacking data will affect the accuracy of determina- 
tion of the phase space distribution function which is the basis of 
the statistical treatment of halo objects. Owing to the small num- 
ber of objects, we will have to assess the influence of the choice of 
halo objects on the accuracy of determination of the {rvl) profile - 
there can be still objects that are not included in the catalogue and 
their addition could change the results. To mimic this situation, we 
will superpose {rvl) profiles prepared for random sub-samples of 
the existing catalogue. 



4.1 The {rvl) profile for compact objects in Galactic halo 

We excluded from the analysis Leol with escape mass > 10^^ M©. 
The high speed of Leo I might be due to a three-body ejection 
mechanism as noted bv lSales et al.l (l2007h . the same remark could 
be applied to some of the fast-receding satellites. Assuming that 
Leo I was gravitationally bound to Galaxy would make the whole 
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Figure 8. Various observables averaged over spherical shells as functions of 
spherical radius. First column for /(e, e) oc ^3/2+1/12 (gi/8 _ ^^^8 
with support 1/12 < u < oo. Second column for /(e,e) = 
^(l - 26^)^ + 144e^ei^ et(l-4e2e)^ and support 1/24 < 
u < 24. 



3.5X10" 




Figure 9. Halo objects in the (r, ^rv'^/G) plane. The values show min- 
imal mass required to gravitationally bind a particular halo object. The 
only object not shown is Leo I for which -^rv^ = 1.03 x IO^^Mq 
(at r = 254kpc). 



work of inferring Galaxy mass from Jeans theory meaningless as 
then Galaxy mass could be decided to be not lower than 10^^ M© 
based on a single measurement, regardless of what value indicated 
the motion of tenths of other objects. As to another object Leo T, 
it cannot be included in preparing the {rv'^) profile owing to the 
fact that it is isolated from other objects on the axis of galactocen- 
tric distance (one cannot calculate dispersion for a single element). 
This fact will not reduce the Galaxy mass estimate due to moderate 
value of escape mass for Leo T which is 1.24 x 10^ \ The remaining 
objects were used to find the {rv'^) profile as described below. 

To determine function (rVr) and see how the choice of objects 
could influence it, we chose 7001 subsamples, each consisting of 
randomly chosen ~ 2/3 of the totality of objects. A particular sub- 
sample was an ordered list of pairs {r,rVr). Each of these lists 
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Figure 10. Distribution of halo objects in space: projection onto the XY 
plane (the Galactic plane) [top], projection onto the XZ plane [middle], and 
projection onto the YZ plane [bottom]. 



was split into sublists of equal length k called fc-boxes. A k-box 
contained objects with similar galactocentric distance. The mean 
values ((r), (rVr)) were calculated for every k-box in a list. The 
k could be chosen small without affecting the reliability of finding 
(rVr) at smaller radii where the concentration of halo objects was 
large - there were hundredths of such /c-tuples at comparable mean 
radii. Averaging over these /c-tuples would give correct result, the 
same as if obtained with a much greater k. Taking small k was jus- 
tified by the fact there was only a few of catalogued objects at large 
radii - too large k would require calculating dispersion for a /c-tuple 
containing halo objects with unacceptably large radial spread. We 
found /c = 6 to be optimal (we used also other /c's to compare with). 
Next, all the lists of averages were joined to form a single output 
list (for k = 6 it consisted of ~ 10^ pairs ((r), {rv^)) - a cloud 
of points on the (r, rv^) plane). Then, the list was sorted according 
to increasing radii and converted to a shorter list with the help of a 
smoothing procedure with a fixed width window (we used 8kpc for 
the width). The data embraced by the window were averaged using 
a Gaussian filter (we assumed the window size to have been 4 times 
greater than the a parameter in the Gaussain function). In effect, we 
obtained the required velocity dispersion profile (rVr) as a function 
of r. It could have occurred that several identical /c-tuples appeared 
in the output list, especially at larger radii, where the number of ob- 
jects was small. Such a /c-tuple was more frequent and effectively 
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it entered with a greater weight into the average. It could have also 
happened that some k- tuples comprised objects with radii very dif- 
ferent form each other, but this situation was rare at lower radii and 
would not affect the final result of the scanning procedure. It was 
more frequent only in the remotest region where the number of ob- 
jects was small. A standard deviation from the mean was used as 
a measure of the uncertainty of determining function {rv^). The 
discussed box-sampling had also the advantage that the effective 
radial window size adapted automatically. For all /c-tuples with the 
mean radii embraced by the scanning window, we obtained some 
effective window width changing smoothly with radius depending 
on the average radial spread of objects in a given k-box. 

The result of the averaging procedure is shown in figure 
Fig.[TT]for several choices of k. It is seen that {rvl) is independent 
of k for small radii and it changes with k for larger radii. Function 
{rVr) is thus reliable at lower radii where the number of objects per 
distance is high and it is much less reliable at the outermost radii 
where only a few objects are present scattered over a very large in- 
terval of distances. This is the nature of the measurement data and 
the situation cannot be improved unless new halo objects are found. 
By comparing with the isotropic model 4(r^;^) / G discussed in sec- 
tion |2.2l we conclude that estimates of asymptotic mass will be sen- 
sitive to the assumptions made about which objects are not gravita- 
tionally bound to Galaxy. With only Leo I and Hercules considered 
as unbound, one can expect the mass to be 1.25 times greater than 
that obtained when, in addition to Leo I and Hercules, also Pal 14 
is assumed unbound. Mass estimates can change also depending on 
k by less than 10%. 



4.2 Galaxy mass in the point mass approximation and a 
systematic way of obtaining the corresponding phase 
space distribution function 

Function {rv^) is a functional of phase space distribution function 
/(e, e). Function / is defined on a triangular support (call it S) as 
shown in figure Fig.[6l but only a fragment of it, an r-dependent 
rectangular region (call it Sr) is used to find {rvl) at a given radius 
r. Function / is not known in advance, it must be found based on 
{rVr) known from observations. This task can be solved in many 
ways. Since / must be nonnegative, it can be expressed in terms 
of an auxiliary function /i(e, e) so that f = > 0. A possible 
idea to find / is to decompose /i in a series of orthogonal functions 
on S. In a first approximation, only several first terms can be taken 
into account. Given a set of expansion coefficients, the correspond- 
ing (rVr) can be calculated by numerical integration over Sr and 
compared with the observed (rv^) profile. A norm of the differ- 
ence between the theoretical and observed profiles (rVr) gives rise 
to some error function of the unknown expansion coefficients. The 
optimal set of the coefficients can be found by minimizing the error 
function. A detailed d escript ion of the minimizing scheme can be 
found in Bratek et"al] (l2012h . In finding / we assumed the radius 
of the internal boundary sphere as small as possible: ra ~ 4kpc, 
and we verified that the mass is not sensitive to this choice as far as 
ra < 24kpc. The radius of the external boundary sphere was set 
equal to the endpoint of the observed (rVr) profile: rt ~ ISOkpc. 

The results of applying the Keplerian ensemble method to the 
measurements of compact objects in Galactic halo (using the min- 
imizing scheme) is shown and described in Fig.fTJ] Assuming that 
Leo I and Hercules are not gravitationally bound to Galaxy (in 
which case the (rVr) profile is prepared discarding the two objects). 
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a) The radial velocity dispersions (r, {rv'^) / (GM)) for halo compact 
objects prepared assuming averaging 6.8, which were the basis to obtain 
the phase space distribution function /(e, e) 
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b) The (r, {rv'^ ) / (GM)) profile [ thin line], a standard deviation from the 
profilQ [dotted line], the model curve for the (r, {rv'^) / (GM)) profile and 
its uncertainty [filled circles], and the prediction for the dispersion profile 
of transversal velocity (r, {rv'^) / (GM)) [empty circles], both derived 
from the found distribution function f(e,e). 



galactocentric distance [kpc] 



galactocentric distance [l<pc] 



c) The resulting expectation value of the symmetrized flattening parameter 
K, of the velocity dispersion ellipsoid: hi((3) = ^ arctan (1 — ^) — 1 for 
a set of orbits crossing spherical shells of various radii. 



galactocentric distance [kpc] 



galactocentric distance [kpc] 



d) The resulting expectation value of eccentricity for a set of orbits cross- 
ing spherical shells of various radii. 

Figure 12. The results of applying the Keplerian ensemble method and the 
minimizing scheme of finding the phase space distribution function to mea- 
surements of Galactic halo compact objects: Leo I and Hercules excluded 
[left: column]; Leo I, Hercules and Pal 14 excluded [right column]. 



we obtain the following estimate of Galaxy mass: 

]\LGai = 2.07 ±0.3 X lO^^M©. 

The mass is slightly lower than necessary to gravitationally bind Pal 
14. However, due to measurement errors of velocity and distance of 
Pal 14 and the uncertainty of McaU Pal 14 is presumably bound 
to Galaxy. If Pal 14 was assumed not bound to Galaxy, the same 
method would lead to Moai — 1.97 x 10^^ M©. By assuming 
that only Leo I is not bound, we would obtain Moai = 2.3 x 
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Figure 11. The result of the k-box averaging discussed in the text. The panels show the profiles (r, {rv'^)/G) in units (kpc, Mq) for halo objects obtained 
using /c-boxes and Gaussian window of width w (symbolically, k.w), respectively, for 5.8 [top row], 6.8 [middle row], 7.8 [bottom row], with the exclusion 
of Leo I [left column], with the exclusion of both Leo I and Hercules [middle column] and with the exclusion of Leo I, Hercules and Pal 14 [right column]. 



10 M© which is less than the escape mass for Hercules of about 

3 X lO^^Mr:). 



5 CONCLUSIONS 

We showed that the motion of halo objects regarded as test bodies 
in the Galactic potential well can be made consistent in the point 
mass approximation with Galaxy mass of about 2.1 x 10^^ M©. To 
achieve this result we gave up the constraint of constant anisotropy 
parameter. In our model the anisotropy of the velocity dispersion 
tensor changes with Galactocentric distance (see, Fig.[T2l). To some 
extent this is associated with variable expectation value for the elip- 
ticity of orbits intersecting concentric spheres of various radii. In 
our method we did not assume the phase space distribution func- 
tion, we focused on finding it based on the radial velocity dispersion 
of compact objects in Galaxy halo. Finding the function from this 
single observable might be non-unique. Our result could be tested 
further. To this end we presented in Fig. [12] the prediction of our 
model for the transverse velocity dispersion. Note, that total mass is 
an asymptotic notion, and dispersion functions are defined locally 
for all radii. Thus, whereas the mass estimate can be reliable, the 
behaviour of observables can be different at smaller radii, in which 
region one should expect a correction from higher multipoles of the 
real gravitational field of an extended mass distribution. 



REFERENCES 

Ackermann M., et al. 201 1, Physical Review Letters, 107, 241302 
Aprile E., at al. 2011, Phys. Rev. Lett, 107, 131302 



Aumer M., Binney J. J., 2009, MNRAS, 397, 1286 
Bahcall J. N., Tremaine S., 1981, ApJ, 244, 805 
Battaner E., Garrido J. L., Membrado M., Florido E., 1992, Na- 
ture, 360, 652 

Belokurov V., Walker M. G., Evans N. W, Faria D. C, Gilmore 
G., Irwin M. J., Koposov S., Mateo M., Olszewski E., Zucker 
D. B., 2008, ApJL, 686, L83 

Belokurov V., Walker M. G., Evans N. W, Gilmore G., Irwin 
M. J., Mateo M., Mayer L., Olszewski E., Bechtold J., Picker- 
ing T., 2009, MNRAS, 397, 1748 

Brainerd J. J., , 2007, The Astrophysics Spectator, [online] 

Bratek L., Sikora S., Jalocha J., Kutschera M., 2012, in prepara- 
tion 

Brown W. R., Geller M. J., Kenyon S. J., Diaferio A., 2010, AJ, 
139, 59 

Clewley L., Warren S. J., Hewett P C, Norris J. E., Evans N. W, 
2004, MNRAS, 352, 285 

Dehnen W, Binney J., 1998, MNRAS, 294, 429 

Dohm-Palmer R. C, Helmi A., Morrison H., Mateo M., Ol- 
szewski E. W, Harding P., Freeman K. C, Norris J., Shectman 
S. A., 2001, ApJL, 555, L37 

Eisenhauer F, at al 2005, ApJ, 628, 246 

Francis C., Anderson E., 2009, New Astronomy, 14, 615 

Geha M., Willman B., Simon J. D., Strigari L. E., Kirby E. N., 
Law D. R., Strader J., 2009, ApJ, 692, 1464 

Geringer-Sameth A., Koushiappas S. M., 2011, Physical Review 
Letters, 107, 241303 

Harris W. E., 1996, AJ, 112, 1487 

Honma M., Sofue Y., 1996, PASJ, 48, L103 

Jalocha J., Bratek L., Kutschera M., Skindzier R, 2010, MNRAS, 



10 



407, 1689 

Jalocha J., Bratek L., P^kala J., Kutschera M., 2011, ArXiv e- 
prints a revised version accepted for publication in MNRAS 

Koch A., Wilkinson M. L, Kleyna J. T., Irwin M., Zucker D. B., 
Belokurov V., Gilmore G. R, Fellhauer M., Evans N. W., 2009, 
ApJ, 690, 453 

Martin N. R, Ibata R. A., Chapman S. C, Irwin M., Lewis G. R, 

2007, MNRAS, 380, 281 
Mateo M. L., 1998, ARA&A, 36, 435 

Sales L. V., Navarro J. R, Abadi M. G., Steinmetz M., 2007, MN- 
RAS, 379, 1475 

Sikora S., Bratek L., Jalocha J., Kutschera M., 2011, ArXiv e- 
prints 

Simon J. D., Geha M., 2007, ApJ, 670, 313 

Sofue Y., Honma M., Omodaka T., 2009, PASJ, 61, 227 

Starkenburg E., Helmi A., Morrison H. L., Harding R, van Wo- 

erden H., Mateo M., Olszewski E. W., Sivarani T., Norris J. E., 

Freeman K. C, Shectman S. A., Dohm-Palmer R. C, Frey L., 

Oravetz D., 2009, ApJ, 698, 567 
Watkins L. L., Evans N. W., An J. H., 2010, MNRAS, 406, 264 
Wilhelm R., Beers T. C, Sommer-Larsen J., Pier J. R., Layden 

A. C, Flynn C, Rossi S., Christensen R R., 1999, AJ, 117, 2329 



